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ABSTRACT 



The Seiico Jo resistance and self -propulsion model test results are 
represented by curve fitting methods which <n functions containing 

Chebyshev polynomials. The accuracy of the curve fitting method is dis- 
cussed, and the range of validity of the functions is presented. 

A computer program for the calculation of EHP and SEP for single- screw- 
merchant ships with hull proportions falling within the range of the 
Series 60 is discussed in detail. Several examples in the use of the pro- 
gram are presented. 

It is concluded that estimates of EHP and SEP can be calculated to 
within a tolerance in the neighborhood of 1$ on the average, and with great 
certainty, to wiuhin a tolerance of less than 3$ for any given speed. 
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I. INTRODUCTION 



1. Foreword 

The determination of the required horsepower necessary to fulfill 
the speed requirements of a new ship design is one of the many problems 
which faces the naval architect in the early stages of ship design. 

He must have some way of knowing what effect various choices of size, 
proportions, and fullness of the ship will have on the power required to 
propel the ship at the speeds he desires. 

In determining this power requirement, the naval architect can draw 
upon his past experience in ship design and the results of specific isolated 
model tests which are published periodically in the various technical 
journals of the profession. 

However, thm results of experiments conducted on families of models, 
in which the different design parameters are varied systematically, are the 
naval architects' most valuable sources of reference in estimating the 
power required for his new design. The result of these tests are usually 
expressed in the form of design charts from which the naval architect, 
using interpolation techniques where necessary, can calculate the power 
requirements for a number of combinations of design parameters, and thereby 
make an intelligent choice of the best combination of parameters to give 
minimum power with due consideration for the other limitations of his design 
conditions. The Series 60 is such a family of models designed to aid the 
naval architect in the estimation of power requirements for single-screw 
merchant ships. 

2. The Series 60 Model Tests 

The research on Series 60 was carried out at the David Taylor Model 
Basin (DTMB) of the United States Navy. Experiments were conducted as part 
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of the Bureau of Ships Fundamental Hydromechanics Research Program during 

the years 1948 to i 960 . The results were published in a series of papers 

presented before the Society of Naval Architects and Marine Engineers, 

and a final report, incorporating all of the information contained in the 

( 1 

interim reports was published in July of 1963 . 

After a survey of current practice in the shipbuilding industry to 
ensure, as far as possible, that the series would cover the normal range 
of proportions of modern ships, proportions were chosen for five parent 
models of the Series 60. These proportions are listed in Table I below. 



TABLE I 

Hull Proportions for Five Parent Models of Series 60 



Model No. 


4210 


4218 


•4221 


4213 


4214 


°b 


.60 


.65 


.70 


• 75 


.80 


l/b 


7.50 


7.25 


7.00 


6.75 


6.50 


b/h 


2.50 


2.50 


2.50 


2.50 


2.50 


A/(l/ioo ) 3 


122.0 


141.4 


163.4 


188.2 


216.5 




6.165 


5.869 


5-593 


5.335 


5.092 


LCB as fo of L 

from 0 

M BP 


1.5 aft 


1.54 aft 


0.55 aft 


1.5 fwd 


2.5 fwd 



For any one block coefficient and LCB position, a total of nine models were 
tested in which the L/B and B/H ratios were varied, for a total of forty-five 
models. To cover the general spread of L/B, B/H, and A/(L/lOO) for 



* 

Numbered references throughout the text are listed on page 44. 
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existing designs and the possible variation in LCB position, a 
adopted as shown by the dotted lines in Figure I. 



FIGURE I * 

Variation of Prop ortions with C K for Series oO Models 



f AFT 




* 



grid 



From reference (l), page II-3 



The pattern for a typical case (C^ = .60) is shown in Figure II. 

FIGURE II * 

Typical Variation of L/B and B/H Ratios for 
a Given Value of 




The midship section of the models had no deadrise in accordance with 

current practice. The relationship between C and C for the Series 6 0 is 

m p 

as follows: 

c = 0.0857 c + .925 

m p 

And, since C. = C x C 
1 b p m 

= 0.0857 C p 2 + .925 C p (1) 

LWL for the Series 6 0 hull forms is related to L by the following 
simple formula : 

LWL = 406.75 x L (2) 

400. 



* From reference (l), page II-3 
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The results of the Series 6 0 resistance tests are presented in the 
fora of easy to use design charts, and in numerical point-by-point fora 
in tables. . 

Values of total resistive coefficient (C ) for forty-five ships of 
length betveen perpendiculars of UCO feet, and varying values of C , L/B, 
and b/H are presented in Thbies B-l through B-45 of reference (l). 

Utilizing the Froude assumption that the total resistance can be divided 
into two parts, the skin friction of an "equivalent plank" and the 
residuary resistance, the values of include a skin friction coefficient 
(C^,) calculated for salt -water at 59°F in accordance with the American 

'ft 

Towing Tank Conference (ATTC) 19^-7 line , and an additional allowance for 
ship correlation (AC^) of + 0.C004- as recommended by the ATTC 19^-7* 

Ifethods for scaling the values of C for an arbitrary ship length '..’ill be 
discussed in Chapter IV. 

Values of propulsive coefficient (PC or EHP/SH?) for forty-five ships 

of length betveen perpendiculars of 6C0 feet, and varying values of C , 

b 

L/B, and B/n are also presented in Tables B-l through B-^5 of reference (l). 
Unlike the values of C , the propulsion data can not be corrected for 
variation in ship length. The choice of a 600 foot ship to illustrate the 
propulsion tests vas made principally because it was considered more rep- 
resentative of modern ships than the hOO foot length chosen for the resistance 
. ** ^ 

presentation. The propulsion data presented vas obtained, using a standard 
propeller diameter equal to 0,70 of the designed lead draft. 

Data for the effect of variation of LCB position on power requirements 
is by no means complete. Corrections to total resistance for variation, of 

ft 

Also known as the Schcenher line. 

rft 

Reference (l), page VII-k. 



- 5 - 



LC3 position for five values of C are presented in 'Tables 49 through 

b 

53 in reference (l), with only four or five different positions of LCB 
noted. 

Although the propulsive data is presented for ships of 6 00 foot 
length, if the designer adopts the lines of Series 60, the position of 
LC3 as used in the parent forms, and a propeller having the standard 
ratio of diameter to draft of 0.7> he can make a very accurate estimate of 
both the EHP and SEP of a ship for any particula.r selection of L, 3, H, 
and A . 

3 . Previous Work in Series 6C Computerization 

The development of a computer program for estimating SEP using the 

( 2 ) 

Series 6 0 results was first attempted by D. Khoushy in 1962. ‘ ' This 
attempt represented the 6 0 by curve fitting methods to some extent 

e.g. wetted surface, but it relied exclusively on interpolative methods 
for Ch which are both cumbersome to program and highly consumptive of 
computer storage space . Furthermore, the large amount of storage space 
required in this type of program made it unsuitable for use as an EHP or 
SEP subroutine in a larger program. It also should be noted that an 
interpolative scheme is not very efficient since the test results are 
not presented in a uniform manner, i.e. at low values of models were 
tested up to a v/Vlw l of 1.10 whereas the higher values of were tested 
up to a V/YLWL of 0.80. This problem was avoided in the Khoushy work by 
discarding a large amount of Series 6 0 data in the low and high ranges of 

v/Vlwl. 

_ _____ _ 

The storage of approximately 6000 numbers is required for the Series 60. 
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) 









Therefore, it was decided that curve fitting vould be the most 
advantageous method of representing the Series 60 model test results 
in a computer program. 
k. The Curve Fitting Approach 

The curve fitting approach for representing the Series 60 results 
is very desirable since once a function is found, which represents the 
given data, one knows the error he can expect in predicting EHP and SHP 
using these derived functions. The fitting of curves through the Series 
GO resistance data was started in early 1966 by T. M. Pitidis-Poutous. 
The basic theory of the curve fitting method employed wa s developed in 
reference (3) in 1965 with the mathematical derivation of ships lines 
being the focus of attention. This curve fitting method, as it applies 
to the representation of the Series 60 resistance data and its extension 
to represent the propulsion results for this work, is described in 
Chapter II. 
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II. DESCRIPTION OF CURVE FITTING METHODS USED FOR THE SERIES 60 



The theory of the curve fitting methods used for representing the 
Series 60 model test results stated herein is largely adopted from reference 

(« 0 . 

Since the chart reading methods of the Series 60 are not suitable for 

computerization, it vas decided to fit curves through the numerical results 

presented for C and PC (or EEP/SHP) in Tables 3-1 through 3-45 of reference 
t 

(!)• 

In the Series 60, both and PC are functions cf four independent 

variables , i.e. V/VEuL, C , L/B, and B/E. Thus, either C or PC can be 

b t 

represented by an equation of the form: 



X T K N 



Y=ima (l) X k K) CL (j^> 

I J * / fc-l Tjr/ J * *> n 71 

and the problem reduces to determining the matrix of coefficients fa n ft J i\ 
from the given data. If the coordinate functions § , , X and n are 

taken as ordinary polynomials, the above expression becomes 



J 



(3) 



I J" K N 



y-iz 



■VZZL) 



i*/ Jr| 

It can be shown that there are certain advantages to be gained if the 

coordinate functions are taken as the Chebyshev polynomials, with no 

additional complications. The independent variables V/T/nWL, C , L/3, 

b 



* See reference ( 5 ). 
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end B/H can be normalized in the range of -1 to 4-1 by the equations given 
in Table II, and the corresponding symbols used for the original, normalized, 
and computer program variables are given in Table III. The equations of 
Table II were derived on the basis of the limits of the Series 60 model 
tests. 

Stipulations on the maximum powers I , J, K> and II, to be used in 

equation ( 3 ) are determined by the amount of data available from the 

Series 60 results. For example, each model was tested for three different 

values of B/H and L/B, therefore I and J must be equal to or less than 3* 

Models were tested for five different values of C , therefore K must be 

b 

equal to or less than 5* Thus, the number of terms in the equation will 
be equal to or le s s than ^5N. The value of h canuuc. be arrived at in the 
same way as I, J, and K, Since each model was tested for twenty to thirty 
different values of v/YLWL, differing from model to model. If we assume, 
for the moment, that N = 10 would be adequate, we see that the determina- 
tion of the coefficients £ would involve the solution of ^50 

linear algebraic equations. This is virtually impossible, both from uhe 
numerical and computer storage points of view. 

These considerations prohibited the use of a direct, parameter fit, 
and indicated the necessity of a step-by-step technique. 

Using the symbols for the normalized variables as defined in Table III, 
equation (3) becomes: 

&n%h)* k (*)n A (r) 

or i i * " 

Y z Y. Y. A fin ( V ) (M 

I J = l 
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table ii 



Equations for Norralizing Independent 
Variables in Series 60 Range 



v = 2.5 - V — - 1.875 

Vlwl 



u = L/3 - 7.0 

1-5 



z = 2.0 B - 6.0 



H 



V_ 
LWL 



= O.kv + 0.75 



= i.5u + 7-0 



B 



B 



H 



= 0.5z + 3.0 



x = 10.0 C - 7.0 
b 



C 

b 



= O.lx +0.7 
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Correspondence of Original, Normalized, and Computer Program Symbols 




> 



> > 




pq 



pq 



n pq 



o 



s 



SI 

& 



o 



s 

o 





H 






O 




rH 


rg 


i — 1 


o 




o 


f 


CO 


f 


CO 


'C* 


CO 




o 




H 


C4 




g3 


*H 


a> 


d 


H 


-p 


•H 


(jjj 


P 


to 


g 


& 


•H 




S 


U 


O 


o 


O 


S 


o 



11 



where : 



* - j 



j( x-. - ) + i j e*. ’) 



and Fg ^ $;(*) IpjM X fe (x) 

The above substitutions reduce the four-dimensional matrix "a" to 
a two-dimensional matrix " A ", which is easier for computations, and 
defines the problem more appropriately for the method to be used. For 
given values of u, z, and x, equation (k) can be written as: 

V P = ^ C en (■*■) ( 5 ) . 

T\ - I 

with L 

C/>y> = ^ ^ ( M > ") > *} (6) 

Equation (5) represents the variation of Y with respect to v for 
the p'kk model, and the IT coefficients (Cp a , n = 1, N) can be evaluated 
from the data in Tables B-l through B- ^5 of reference (l), using the 
least-squares criterion: 

Z - Z (><:)} 

by differentiating partially with respect to the unknown coefficients. 

All models w r ere fitted in this fashion, and the resulting equation 
for each model was evaluated at intermediate values in order to ensure 
that the polynomial had no undesired oscillations. In fitting the curves 
for C , it was found that a reasonably good fit was obtained with I! = 10 

L 

(corresponding to a 9th. order polynomial) and since this is rather high, 
the possibility of introducing a suitable weighting function was explored 
Therefore, the fittings were repeated with the following two alternative 
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weighting functions assigned to C : 

t 



Y - C<c*x/o 3 ) / vvrr^iy 



The last alternative was found to' give the best fit with N = 10 for most 

models. Going from N = 10 to N = 11, the maximum error experienced in 

C x 10^ only decreased from .025 to .024 for Table B-l of reference (l). 

Having obtained matrix (cl corresponding to equation (5) with 

p = '45 (number of models), matrix Oml was obtained by exact fit, i;e. 

L = 45- It should be mentioned here that a certain amount of regression 

could also be tried on this step of fitting, and this was actually 

attempted, but without success. The reason is that, although it might 

be possible to obtain a good fit for C or PC at a constant value of 

t 

V/lfLWL with L less than 45, it does not follow that such a fitting could 
be applied successfully to the coefficients of the C - v/VLWL equation 
(i.e. C^). Thus, direct 4-parameter fitting of the data could have 
produced an equation with less error and less terms, if it were possible 
numerically. 

The values of PC were fitted in the same manner described above, but 
no weighting function was used. 

The final equations obtained are as follows : 



C* x /o 3 



10 4f 



(V/Vz^n.) 



(7) 



/o 45" 



pc= fr II 



*--■1 



(8) 



-13- 



where 



F, - T:-, («) Tj„ ( 5 ) T ft „ (*) 

Tp — p ^ Chebyshev Polynomial 

The elements of the matrices IA land |B \ are read in as data 

1 nl* 1 nl J 

in the main program, and are listed in Appendix A. These matrices form a 
compact computational package of $00 numbers for the calculation of C^_ 

(400 ft. ship) and PC (600 ft. ship). 

No attempt vas made at fitting curves through the data for correction 
to total resistance due to variation of LCB position. The inadequacy 
of the data presented, as veil as the erratic nature in which it behaves 
with v/VlWL, precluded the use of the curve fitting methods described above. 
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III. RESULTS OF CURVE FITTING METHODS USED AND RANGE OF VALIDITY 



The results of the curve fitting methods used for C and PC were 
first tested by reproducing the date contained in Thbles B-l through B-k5, 
i.e. for the forty-five models of the Series 60. The results of this 
test are presented in Table TV, vhere the maximum error, in percent, 
experienced in C and PC is noted for each of the forty-five models. 

w 

The maximum error experienced in C was 2.7^0$ for the parent model 

4218. The maximum error experienced in PC was 2.030$ for model 4276. 

Although these figures appear prohibitive at first glance, it should be 

emphasized that these are the maximum errors experienced for any given 

model, and they occur at only one value of v/ifLVL. The errors at other 

values of V/YLWL for any given model are well under 1$ in the great 

majority of cases. In fact, as can be seen from Table TV, the maximum 

errors experienced in C and PC for any given model are well below 1$ 

in most cases. To illustrate this point further, plots of C versus 

v/VlWL are presented in Figures III through VII, and plots of PC 

(or EHP/SHP) versus v/VLWL are presented in Figures VIII through XII, 

for the five parent models of the Series 60. In these plots, the solid 

line represents the value of C and PC as calculated from equations (j) 

and (8) respectively, and the points represent the value at that v/VLWL 

presented in the tables of reference (l), which are in error enough to 

be visually discernable on the scale used in making the plots. It is 

further noted that Figure TV represents the worst case experienced in 

fitting either C or PC. 

t 

As a second test, equations (7) and (8) were evaluated by generating 
125 interpolate and extrapolate models, especially concentrating on 
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TABLE IV 



I-kximum ^ Errors Experienced in. and PC for the 45 









Model 


s? of 


the Series 60 




Series 60 


Model 








Max. <p Error 


I&x °!o Error 


Table No. 


No. 


C b 


Ul 


b/h 


in C+. 


in PC 


B-l 


4240 


.60 


6.5 


2.5 


.578 


.853 


B-2* 


4210 


.60 


7.5 


2-5 


1.230 


.648 


B-3 


4267 


.65 


8.25 


2.5 


1.515 


.672 


B*4 


4264 


.65 


6.25 


2.5 


.340 


.616 


B-5* 


4213 


.65 


7.25 


2.5 


2.740 


.298 


B-6 


4243 


.60 


8.5 


2.5 


.701 


.696 


B-7 


4244 


.70 


6.0 


2.5 


• 531 


. .248 


B-S* 


4221 


.70 


7-0 


2.5 


.370 


.402 


B-9 


4247 


.70 


8.0 


2.5 


.232 


.552 


B-10 


4268 


.75 


5-75 


2.5 


.819 


.366 


B-ll* 


4213 


• 75 


6.75 


2.5 


.824 


.697 


B-12 


4271 


.75 


7-75 


2.5 


1.045 


.682 


B-13 


4243 


.30 


5-5 


2.5 


.191 


.372 


B-l4* 


4214 


.80 


6.5 


2.5 


.381 


.134 


B-15 


4215 


.00 


n s 


2.5 


.243 


.394 


B-l6 


4252 


.60 


6.5 


3.0 


.885 


.437 


B-17 


4255 


.60 


7.5 


3.0 


.714 


.393 


B-l3 


4254 


.60 


8.5 


3-0 


.477 


1.261 


B-19 


4272 


.65 


6.25 


3.0 


.602 


.529 


B-20 


4275 


.65 


7.25 


3.0 


.405 


.446 


B-21 


4274 


.65 


8.25 


3.0 


.600 


.516 


B-22 


4256 


.70 


6.0 


3.0 


• 175 


• .277 


B-23 


4259 


.70 


7.0 


3.0 


• 134 


.590 


B-24 


4253 


.70 


8.0 


3.0 


.278 


.573 


B-25 


4276 


• 75 


5-75 


3-0 


.523 


2.030 


B-26 


4279 


.75 


6.75 


3-0 


1-331 


1.430 


B-27 


4273 


.•75 


7-75 


3.0 


.923 


.591 


B-23 


4260 


.30 


5-5 


3.0 


.248 


.517 


B-29 


4263 


.80 


6.5 


3.0 


.166 


.129 


B-30 


4262 


.80 


7-5 


3.0 


.082 


.284 


B-31 


4241 


.60 


6.5 


3-5 


.541 


.449 


B-32 


4253 


.60 


7-5 


3-5 


.4ll 


.878 


B-33 


4242 


.60 


8.5 


3-5 


.459 


.906 


B-34 


4265 


.65 


6.25 


3-5 


.634 


•537 


B-35 


4273 


.65 


7.25 


3-5 


.418 


.349 


B-36 


4266 


• 65 


8.25 


3-5 


.493 


.664 


B-37 


4245 


• 70 


6.0 


3-5 


.203 


.443 


B-33 


4257 


.70 


7-0 


3-5 


.435 


.452 


B-39 


4246 


• 70 


8.0 


3-5 


.312 


.544 


B-40 


4269 


• 75 


5-75 


3-5 


..223 


.598 


b-4i 


4277 


• 75 


6.75 


3-5 


.506 


.581 


B-42 


4270 


• 75 


•7-75 


3-5 


.919 


1.275 


B-43 


4249 


.80 


5-5 


3-5 


.192 


.542 


B-44 


4261 


.80 


6.5 


3-5 


.191 


• 577 


B-45 


4250 


.80 


7-5 


3-5 


.118 


.303 



* Denotes parent model 
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intermediate values of the fiveCb's used for the Series 6 0 models. This 

vas done in order to determine the range of validity of these equations 

with respect to their four independent variables. It should be emphasized 

that the ranges of validity determined below are not just a result of 

the inadequacies of the curve fitting methods used; they are very much a 

result of the limits of the Series 60 model tests themselves and the 

range of data made available in reference (l). 

The data presented in reference (l) is based on the use of three values 

of B/H, i.e. 2.5, 3*0, and 3*5* Since equations (7) and (8) rely on a 

second order interpolation scheme in B/h , which is the method suggested in 

reference (l) when using the design charts, there should be no discrepancies 

caused by this variable. It is also expected that a limited amount of 

/ * 

extrapolation is possible in B/H. 

Similarly, there should be no discrepancies due to L/B since equations 

(7) and (8) rely on a second order interpolation scheme for this variable, 

which is also the method used by BTM3 in the preparation of the Series 60 

design charts. However, as can be readily seen from Figure I, there are 

limits (imposed by the proportions of the Series 60 models) on the range 

of interpolation possible in L/B. The range of interpolation in L/B is a 

function of C and can be expressed as follows: 
b 

(L / B) min = 5-5 + 8 -° ~ 10 C b (9) 

2 

( L /3)max = 7*5 + 8.0 - 10 C b (10) 

2 

- • — 

See example of the ship SCHUYLER OTIS BLAND in Chapter V. 
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A small amount of extrapolation beyond these limits should not cause 
trouble . 

In producing the Series 60 design charts, a third order interpolation 

in C vas used by DTMB. In equations ( 7 ) and (8) a fourth order inter- 

polative scheme for C is used^vhich could produce slight discrepancies 

in interpolation and extrapolation. It is recommended that the program 

be used only in the range of C = 0.60 to 0.80. 

b 

The Series cO model tests impose definite limits on the range of 

V/y'LWL permissible. For example, in testing equations ( 7 ) and (8) vith 

the 125 models mentioned above, it vas found that vith = 0.775> B/H = 

and L/B = 7*0, the values of C and PC given by equations ( 7 ) and (8) 

t 

dropped off sharply at v/Vlwl = .925 and e.c.. v.^nt negative at V/VLWL = .950* 
Ihis is understandable biu^e at these values of V/VLWL ve are outside the 
range of results given by the Series 60 tests. Therefore, to avoid 
occurrences of this nature, a "trap" has been built into the EHP-SKP program 
described in Chapter IY, vhich puts limits on the value of V/VLWL a Hove d } 
depending upon the value of C . These limits are presented in Table V. 



TABLE V 

Limits of v/lClWL Alloyed by Program 



Range of C 

b 


(v/Vlwl) . 

N 7 / min 


(v/Vlwl) 

max 


.600 < C < .612 

b ~ 


.45 


1.10 


.612 < c b < .625 


.45 


.90 


.625 < a < .725 


. 4 o 


.90 


.725 < c b < .775 


• 35 


.875 


•775 < C b £ .800 


• 35 


CD 

O 

O 
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Within the ranges of independent variables stipulated above, one can 
expect values of C and FC from equations (7) and ( 8 ) to agree with the 
published results of Series 60 to within Yp, with consistent reliability . 
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IV. DESCRIPTION OF THE 2EP-3H? COMPUTER PROGRAM FOR SERIES 6 0 



The program for calculating EHP and SHP using the Series 6 0 model 
test results was written in the FORTRAN IV language. A flow chart for 
the program, a listing of the FORTRAN statements, and the 900 elements 

of the matrices [A J and IB ,1 mentioned in Chapter II are included in 

L nl J 1 nl J 

Appendix A. The following is a general description of the methods of 
calculation used. Detailed instructions for the use of the program are 
contained in Appendix C. 

1. Inputs 

First the elements of matrices (A , "\ and |B .1 are read in as data. 

1 nl J 1 ruJ 

The number of ships to be calculated (nJGBS) and an identifying alphameric 
statement is then read. It takes three cards to describe each ship to 
be calculated. The first card describes the input mode (iNOPT) for the 
ship to be calculated. The operator is afforded an option of one of two 
possible input modes for each ship, 
a. Mode I 

Entry is made with L, B, H, C^, ^ max > ^min; ^increment > Position 

of LC3 from ^ , S, -f , and kinematic viscosity, 

b. Mode II 

Entry is cade with L, B/H, A, C p , V^, V min , V increK:ent , AC f , Position 
of LC3 from ^ , S, and kinematic viscosity. 

2 . Calculation of needed parameters 

Using the input information, the program then calculates the remaining 
parameters required for executing the remainder of the program. For example, 
if the entry is made using Mode I, the program calculates A, b/e, C p using 
equation (l), the Series 60 equivalent LWL using equation (2), the volume of 
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displacement ^7 , midship coefficient C , and L/B. 

3 . Normalizing hull proportion parameters 

The program calculates the normalized hull proportion parameters 

for L/B, B/H, and C using the eauations in Table II. 

b 

4. Trap to make sure v/YLWL is within range of model test 

The program calculates v/YLWL maximum and minimum and ensures that 
these values are vithin the range of validity given in Table V. If they 
are not, the program alters V rrn - n and V way accordingly, until they are 
vithin the range. 

5. Calculation of number of increments in V to be calculated 

Using V Lrax; v Jn ^ n , and Vj_ ncrerrier) t^ the number of different speeds (m) 
at which EHP and SH? are desired is calculated. There is a maximum of 
100 different speeds allowed for each ship. If this number exceeds 100, 
the program keeps doubling ^increment until "m" is 100 or less. If the 
operator desires EKP-SHP information for only one speed, m is set at one. 

6. Calculation of V's, v/YLWL's and normalizing V/TTLWL 

The V's and v/VLWL' s are calculated based on the Vm n , Vj_ ncrerien ^ 

LWL, and m calculated above. The v/VLWL's are then normalized using the 
equation in Table II. If = V E i n and Vj_ ncreir . 3n t =0., only one speed 

is used. 

7. • Calculation of C„_ for a 400 foot ship and PC for a 600 foot ship 

1 ' 1 ' ' ' ' " ' L> “ . . r - - - ~’ 1 ~ ~ 

The C^_ for a ship of L = 4C0 ft. (C.^) and PC for a ship of L = 600 
ft. are calculated for each v/VLWL -using the normalized parameters determined 
above, and equations ( 7 ) and (8) respectively. Subroutine CEF is used to 
* See. program listing in Appendix A. 
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evaluate the Chebyshev Polynomials for these equations. 



8 . Calculation of C for a 400 foot ship and C <. for a ship of length L 

The Reynolds Numbers (Re) for the L = 400 ft. ship (Re„ ) and for the 

M 

ship of length L (Re^,) are calculated for each v/YLWL by using equations 
(11) and (12). 



406.75 x v s 

1.2317 x 10-5 



( 11 ) 



Re, 



LWL x v. 



( 12 ) 



The C ' s for the 400 ft. ship (C^ ) and ship of length L (Cp_) are 
f 1 M 3 

calculated using the ATTC 19^7 correlation line given in equation (13). 



-•5 



= 4.132 log 1Q (Re x C f ) 



(13) 



In the calculation of the C^'s, it is necessary for the program to use an 
iteration process, since appears on both sides of equation (13). The 
initial estimate for this process is made by using the equation for the 
International Towing Tank Conference (ITTC) line, equation (l4). 



c; = o>Q7 5 _ (lU) 

(log 1() Re - 2 T 

Equation (13) is recycled for 50 times, or until the C^ values converge to 
within a tolerance of 5«0 x 10 ” ^ . Although this should not occur, if 
formula (13) does not converge, the program prints the error signal 

* Rote that LWL = 4o6.75 for a ship of 400 ft. length between perpendiculars, 
and that kinematic viscosity is for salt water at 59"R- 
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1 









"SCHOEIiHER NO CONVERGE" with the speed at which the trouble was experienced. 



9 . Calculation of wetted surface (S) if not input 

Ideally, the wetted surface (s), in square feet, is fed into the 
computer as input. However, if the wetted surface is not known at that 
time, the program will calculate an approximate wetted surface using 
equation ( 15 ) ♦ 



s = 2.654 (L/B)' 32 ' (b/h)’ 2Cj (c,)- 0009 (\7)' 



.666 



( 15 ) 



This equation for approximating wetted surface for the Series 60 was 
developed by D. Khoushy in reference (2). This formula was tested exten- 
sively and comnared to the contours for wetted surface calculation. Figures 
B-124 through 3-126, in reference (l). It was found that *>is formula 
adequately represents the wetted surface contours of the Series 60, with 
the maximum error rarely exceeding 1.2$. Since wetted surface is a factor 
in calculating resistance, EHP and SHP, if wetted surface is not an input 
parameter to the program, one can expect errors in EH? and SHP due to the 
use of equation ( 15 ), in addition to the errors to be expected due to the 
use of equations ( 7 ) and (8) as discussed in Chapter III. 



10 . 



Calculation of C for the ship of length L, R, EHP, and SHP 

The C for the ship of length L, (C t ^), is calculated using equation 



(16). 



C + = C + - C- - 0.0004 + C_ + AC 
t S Hi f M f S f 



(16) 



Then the remaining calculations are trivial, namely: 
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R = 1/2 f S v 2 C 

s t S 

EHP = R v 

s 

550 

SHP = EHP 
PC 

where PC is for a ship of L = 600 ft. 

11. Calculation of LCB position from $ as a percentage of L 

The LCB position from ^ as a percentage of length between perpen- 
diculars is calculated and presented in the output. It should be remembered 
here that corrections to R, EHP, and SHP for variation in LCB position from 
those of the parent model hull forms are not built into the program. Graphs 
for assistance in making these corrections, based ou To 1 ' !!'''• U 9 through 53 
of reference (l), are included in Appendix B. The LCB position from ££ 
as a percentage of length between perpendiculars is calculated for entry 
to these graphs. 

12 . Output 

The output for each ship calculated includes C^, L/B, B/H, C , Cm, 

length between perpendiculars (LBP) in feet, LWL in feet, 3 in feet, H in 

feet, displacement (DISP) in tons, volume of displacement (VOL) in cubic 

feet, wetted surface (WSURF) in square feet, water density 'f in 

lb. sec , kinematic viscosity (GMJ) in f tr , LCB position as a percentage 
ft. 4 sec 

of LBP, and a tabulation of V, V/VlI/L, R, EHP, PC, and SHP for each speed 
requested in the input. A statement is printed under this tabulation which 
reminds the user that the values are not corrected for variation in LCB 
position as explained above. 
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V. EXAMPLES IN USE OF THE EHP-SHP COMPUTER PROGRAM FOR KNOWN SHIPS 



The final EHP-SHP computer .program was tested with the parameters 

of six ships whose EHP-SHP characteristics are known. As a first test, 

the Series 60 equivalent of the ship SCHUYLER OTIS BLAND was chosen, 

since this was the ship used in Appendix D of reference (l) in illustrating 

a numerical example of the use of the Series 60 charts. It was also 

chosen since it illustrates the ability of the program to interpolate for 

intermdeiate values of C (0.651) and l/b (6.82) and to extrapolate 

b 

outside the range of the Series 60 in B/H (2.444). 

As a second test, the five parent models of the Series 60, with 
L = 600 feet, were chosen since complete information on their EHP and SEP 
behavior is presented in Tables 27 through 31 in reference (l), and they 
cover a variety of hull proportion combinations. 

The results of these tests are discussed below. 

1. Test using the ship SCHUYLER OTIS BLAND 

As a first test on this ship, input Node I was used with the wetted 
surface area of 39994.0 square feet as input. As a result, the computer 
output. Table VI, for EHP was a measure of the accuracy of equation ( 7 ). 

In comparing the computer program's EHP results with those resulting from 
the hand calculations carried out in Table D-3 of reference (l), it was 
found that the maximum error experienced in EHP wa s + 2.70$ at v/VLWL = . 85 , 
while the average error experienced was + 1.19$ • Although the average error 
of 1.19$ is fairly gratifying, the maximum error certainly is not. However, 
referring to Table T7 again we see that in using the hull proportions of the 
SCHUYLER OTIS BLAND we place ourselves in the neighborhood of parent model 
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Computer Output SCHUYLER OTIS BIAM), Wetted Surface as Input 
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Computer Output SCHUYLER OTIS BLAND, No Wetted Surface as Input 
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k2ld where the maximum error of 2-7^$ was experienced in the test of 
equation (7)« 

The SCHUYLER OTIS BLAKD ras also tested without wetted surface as 
input. The results of this run are presented in Table VII. We can see 
that the wetted surface calculated by equation (15), 39^51*5^ square feet, 
is in error by -.85850. The error in the calculation of wetted surface area 
tends .to reduce the error in EHP in this instance. 

One must keep in mind that we are comparing -computer output with 
values derived from hand calculations, in which design charts were read with 
an unknown accuracy. 

2. Test using the five parent models of Series 60 with L = 600 feet 

The results of the computer run with the five parent models of 
Series 60 with L = 600 feet are presented in Tables VUE through XII. In 
these cases, the wetted surface area of the ships was calculated by equation 
(15), since it was not known. The maximum and average percent errors in 
EHP and 3K? experienced in this test are presented in Table XIII. Although 
at some points the errors exceeded 2$, in no case did the average error 
exceed 1$. These errors are probably due mainly to the innacuracy of the 
equation used for wetted surface area. 
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VI. CONCLUSIONS AID RECOiSdSNDATIONS 

Based upon the results of tests run to date, one can conclude that 
the EHP-SHP computer program presented herein can be relied upon to 
estimate the EHP-SHP requirements of single-screw merchant ships, falling 
within the range of hull proportions of the Series 60, to a tolerance in the 
neighborhood of ip on the average, and with great certainty, to within a 
tolerance of less than 3t> for any given speed. 

The weak link in the EHP-SHP program appears to be the calculation 
of wetted surface area. It is recommended that the representation of the 
Series 60 wetted surface contours be further investigated, to determine 
whether equation (15) can be modified and improved. 

It is further recommended that the elements of matrices [A 1 and 

nl' 

(vl be punched in the form of FORTRAN data statements, and included in the 
body of the rain program. This will greatly facilitate its conversion to 
a subroutine for inclusion in a larger program if it were so desired. 

It is also recommended that the program be extensively tested with 
existing ships, to further determine its accuracy. 
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1, Flow Chart for FHP-SHP Program 
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READ : 
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rwEAti : 
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Calculate : 

I>‘* L , V , Cj« , L/3 

Normalize : 

L/3, B/H, C b 




— ^SAX/VdFl : VLMAa ^- 
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Calculate : 
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2. Listing of EHP-SHP program and S/R CEF 



C PROGRAM TOC-.LC EMP,SHP USING SERILS 60 MODEL TcSI RESULTS 

01 4 FUSION A ( 10,45) » 3 ( 1C, ^ 5 ) 

L) I M E M S I ON VS( i DO ) , VL( 100 ) ,V ( IOC ) ,AVL ( 13 ) , 3VL ( 1C ) , PHI ( 10 ) 

D I •! E M S ItJ'-j PS I C L' 1 ) , CM I ( 1" ) , v/ P H I { 10 ) , C M( L TO ) , PC M ( 1 00 ) , C F S ( 1 3 0 ) 

l) I 'MENS I ON CM( LOO) ,CTS( 100) ,R( 100) ,EH?( 100) ,SHP( 100) ,QJ( 20) 

L EUR -I AT ( 5F16. 7) 

2 FORMAT ( I"-t LlA't) 

3 FORMAT ( 7 FLO. 2 ) 

4 FUR 'I AT (4F1D. 6,610.6) 

5 FURMaT ( 26H SCHUcVIER N) CUNVtRGE. VS=»F6.3) 

6 FOR MAT ( 4 riO C> — j Fb.3» 2X» 5h L/5=,F6.3,2X,5H t>/H=,F6.3»2X»4H C- = i 

1F6. ) , 2 X » 4H C M = , F o . 3 ) 

7 FOR-UT ( 5H L3P=,Fd.2, 3H FT,2X,5H LwL=,F8.2»3H FT,2X,3H 8= , 

1 F 3 . 2 » 3 H El , 2 X , 3H H=,F6.2,3H FT) 

6 C URMAT ( 6H U I SP = , F 10 . 2 , 5rt TjNS,2X,5H VOL= ,”F 12 . 2 » 6H FT**3‘,2X,' 

17H WSJK F=, F lo. 2 , cH FT**2) 

9 FORMAT ( IH0.2X, 7MV '( X T S ) , 3X , 7H V /R T L /, L , 6 X , 7HR ( LBS ) » 5X » 6HEHP £ H ~ 1 • 
L4X, 2HPC, 7X, 3HSHP (HP)) 

10 FoR MA<- ( F 9 . 3 , F 1 0 . *+ , F 1 4 . 2 , F 1 3 . 2 » F b . 4 , F 1 3 . 2 ) 

500 FORMAT ( TH LC2 IS,FS.3,33I PERCENT CJ«= L39 AFT OF AMIDSHIPS.) 

' 50 L FORMAT’ ( 2 1 H L CO IS AT AMIDSHIPS. ) 

502 FORMAT 17H LC6 IS,F6.3,33H PERCENT Or L3P FWD OF AMIDSHIPS.) 

503 FUR MAT (14) 

5C4 FOR AAT ( 1H1 , 14, L9A4) . 

505 FORMAT ( 1CH DELTA CF= , F 10 . 7 , 2/ , 5h RHO = , F 1 0 . 7 , 1 6 H L 3*S EC**2/ F T *= - 
1 2 X , 5H GNU=, Elb.9, 1CH FT**2/$EC) 

5C6 FORMAT ( 74H0 AMOVE VALUES ARE ~NO T CO'RRECTED FOR ANY' VAR I AT > I 

1 JF LCl POSITION ***) 

C " RtAO COEFF OF POL FOR CTM AND PCM 
00 20) N= 1 , ID 

READ (5,1) ( A ( N » L ) , t. = 1 , 1 5 ) 

2 CO CO'IT I DUE 

DO 201 M=l , 10 

Ac AD (5,1) ( o ( N » L ) , L = 1 , 4 5 ) 

231 CONTINUE 

C READ N.J'IDER oF SHIPS AND ID I.nFO 

"READ (5,2) NJ03S, (00(L) ,L=1, 1 n 
DU LOUD N J IM = I , N J U 3 S 

C READ I ii PU T A )D CaLC NEEDED P^R.^ JETERS FOR SHIP 

READ (5,503) UUPT 
IF (IrJUPT) SOC, 9 00, VC I 

90 C READ ( t , i ) XLbP, 3H, OE L , C P , V U l J , V M A v , V I NC 
C ii= . 0 :>5 7* ( C P ) *"2 + .925 4- C P 
H = S0RT ( ( 35.'.-uEL ) / ( C 3 * X L .■> P * o H ) ) 

6 M = u H * N 
GU TU 933 

Sul R cA D ( 5 , 3 ) a L » P , o *1 , H , C , V . I , , V *U\ \ , v t j ^ 

DEL = C n*XL 3 P* * rl / 35 . 
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0 375) 



G J 3 



C 



O . 



c 



60 0 



£C 1 
c 2 2 



oC 4 

1 1 
12 
13 
l*t 
C 

1 5 



16 



dH= 3 0/ H 

C ? = ( - . 9 2 5 + SOPH .925**2 = 4. 
READ 15,h ) JCF * \LC.-i , GS'JxF , nIJ, 
X L .. L = 4C 5 . 7 o * X L 3 J / 4 0 f . 

700=35. *0EL - - - 

C -l=C5/CP 

3L=XL;i?/J.i . 

;;j \ -iAL i z e pa.-uoe teas 
U = (- 7.- tiL 1/1.5 

Z~~o . 0 f 2 . ~ *BH 
• X =- 7 . 0 * 10.3*03 
T PA P TO *:AKi= SJ3E VL 



.375*0 3 ) ) / ( 2 . 



■; J 



IF ( C 6 . o r . C 
IF (Cd.GT.2 
IF ( C3.GT.C 
IF (C6.GT.C 
IF (co.or.o 

V LO I N= . 4 5 

V L A X = 1 . 1 
GO 10 11 
VLM I •*= .45 
V L‘-i AX=. vO 
oU To 11 

;l i . = .40 
VLOA.< = .90 
30 TU -1 1 
VL.'. r:= . 3 5 

V L 3i X = . d 7 5 
•oU IU il 

7 L H \ = . 3 5 
7 L i A a = . 5 ^ 



A !\ 0 . O 

AX 

6 2 5 . A N 3 . C 0 
725 .AMO. C 3 
A M 0 » C 3 



IS 

. L 



6 12 . A M 0 . 0 2 



n : 



X IT H I 

.0.6121- G 
. E . 0 : . o 2 5 ) 

. E .0 . 7 25 ) 



RANGE OF 7.00 EL 
TO 300 

TO 



TEST RESULTS 



.Lc 
. LE 



. 775 ) 

.3 00 ) 



uJ 

GO 
G o 



TO 

ro 

ro 



601 
602 
30 3 
604 



IF (V-U 1/IaLaL**.5)-VLMIM) 12, 13, 13 

V .-i I M = 7 L . , I O* ( A L ,i L* * . 5 ) 

IF ( V L ■! A X- V -1 A X / ( X L w L * * . 5 ) ) 1 4 » 15,15 
7 0 A a =. 7 L .1 A X * C X L . > L * * . 5 ) 

CALC OF 4. IF 7 1 .0=0 ., -1 = 1 T.-iXUJGriCuT PHUGAAO 
IF (VI :C . uT . 0 . ) GO TO IS 
0 = 1 

7 -110 = 7 -I AX 
GO TO 20 

7 1= I 7 .-! A X- 7:\ I 11/7 INC + 0.5 
•1 = 7 A 

. 1=0 + 1 

1 A\ E SJO.T -1 MOT GREAT EG THAI l ' 0 
IF lO.LE.lOO) GO Ij 20 

VI . C = 2 . ' - : V I -,C 
GO TO Id 

CALC •- H SHIPS 7, '7 OVck or L L , .u « IAL I 2c 7 JVcR R T LwL 
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-•>0- a = -'- a- 

0 r .f jj [ = L , N 

... vs( i ) = + ?*vuc - ... 

VL( I ) = V S ( I ) / ( XL*L**.5) 

— V 1 I )=2.5*v'L( I )-l.a/5 
?=?+!. 0 
30 CONTINUE 

CALC OF C T ’ 5 A ;D j CS F-)K 400 FT UDtL 

-CALL C£F( 5 ,X,CHI ) - . ... . 

CALL CEP ( 3 » U » PH I ) 

CALL CEF { i, Z , PSI ) 

•JO 32 N=i,iC 

A V L l :M ) =0.0 .... . 

3 VL ( 0) = D.C 

.32 CONTINUE . 

00 3 5 K = 1 , 5 

00 35 1=1,3 

DO 3 5 J = 1 » 3 

L = J + 3*1 I- l J + K-l ) 

JO 35 N=I,10 

AVL l i\i l =AVL t ;w +A ( r: : L i * P n I ( J )*.->$ I (I) *CH 1(2.) 

3 V L ( N ) = 0 V L ( N ) + -> ( N , L ) * P H I ( J ) ^ . i i / ■ 0 ' J ^ 1 * ) 

35 CJsTiNJE ‘ 

Dj 40 I -i = 1 » M 

CALL Ci£r( 1C,V( IN) , VPHI ) 

C r-t( IN) =0.C 

- P C.H ( I.-!) =0.0 

00 36 .'4=1 , 10 

.c ra ( i .-I ) =CT;n i i m ) + avl ( t\ ) * vph n n ) - 

36 ?CM< I M ) =PCK ( IN) +oVL(N)*VPHl ( I) 

C Til ( I A ) = C Tri ( I “0 *V L ( I N ) * V L ( I M).*G . OC i 

40 CONTINUE 

CALC OF CF'S FJ.< 400 FT MODEL AND 002 SHI? LENGTH 

00 50 1=1, M 

- 2E.4= ( VS ( I ) * 1 . o 3 o 0 £ 4 C 6 « 75 ) /( 1.2ai /e-05) 

RES = ( V S ( I )*i .o3iV*/LUL ) / t GNU) 

~r.I( I ) = J . 0 7 5 / ( ( ( ALOG ( ?.£ A ) / 2 . 30 25 35 1 ) -2 . 0 ) ** 

CFS( I ) =0 .075/ ( ( ( A L 0 G ( N E S ) / 2 - 3 0 2 5 c 5 1 ) - 2 • 0 ) * * 

DO 43 K. = 1 , 5 C . _ . 

DJ'1‘iy = CF A( I ) 

3 0 G 0 S = C r S ( 1 ) 

CF 1( I ) = ( 2.3025651/ (4. 132*ALUG( <c.-t*CM( I)) ) )**2 

CFSl I ) = ( 2.30 25 551/ (4. 132*AL JG ( *<= $ *CFS( I ) ) ) ) **2 

1 F ( AES ( JvJi'-i.-l Y-CF •)( I ) ) - 5 . ' E-C 7 ) *t 1 , 4 1 , 4 5 

41 I F( AJS 1 j'JGO S-CF S ( I) ) -5 . .£-07)5^,50,45 . 

4 5 C 1 J 4 T 1 NO E 

■TU Tc l o,5) VS ( I ) 

5 0 C Li -i T I i 0 c 
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r\j rvj 



C CALC -OF •.'I S u R F I 1- NOT GIVE n 

If- ( wSURF . GT .0. ) GO TO 52 

WSJ R r = 2 . 6 54 * l JL ** . 3 2 9 ) * ( uU* * . 20 o ) * ( 1 . /CO** . 0 609 ) * ( VCJL ** . 666 ) 
C CALC OF OUTPUT 

62 00 60 I = L.o - ... - - . 

erst I ) = C T N ( I ) -CFO ( I 1-0.0 0 0 A +C F S ( D+OCF 
. . R ( I 1=0. 5* RriC* AS JRF * ( ( V S ( I ) * 1 . n do 9 ) **2 l*CTS( I ) — ._ 

ErlP m=(R( I) -VS ( I 1 - 1 .6860 1 /530 . 

SHPC D = tHP( I ) / P C ' 1 ( l ) 

60 C JO T [ imUE 

... --WRITE (6,506) UJUtSS, ( 0 0 ( L ) , L = l » 19 ) _ ...-. 

n R I T - (6,6) C o , 6 L , 3 H , C ? , C M 

wRI TE (6,7) XL3P , XL..L , 3 -1 , H ... 

WRITE (6,0) DEL, VOL, W SURF 

WRITc 1 6 » 5 C 3 ) UC F , R HU * GN ’J - . _ - - 

C CALC POSIT OF LC3 AS PERCENT LbP 

• PLC d=A dS ( XLCt-lOC . / XL u P ) . . 

r F ( XLCB ) 30^ , 301 , 302 

300 WRITE (6,500) PLUG - — 

go T- 

30 l WRITE ( 0,501) -- _ _ 

uL) TO 6 w 

302 .-.RITE ( o , 50 2 ) PLC3 ... .... ... 

8C WRITE ( o » 9 ) 

- OU 70 1 = 1,0 

WRITE (6,10) VS ( I ) , VL ( I ) , h( I ) , E.iP( I ) , PC -1 ( I ) , SUP ( I ) 

70 CGwT I.'iJE . - . - - 

.-.RITE ( o,506) 

— LOGO CJNTI'iJE _ .. _ 

S TUP 



S *35 ROUT I IE Ccr( I J,U, OEV) 

01 -1ENSI JM C c V ( 10) 

CEV l 11 = 1 . 

CEV ( 2 ) =J 
00 l 1 = 3,10 

1 CEV ( I ) = 2 . - 0 * C E V ( 1-1 ) - C E V ( 1-2) 

RETURN _ _. .... . .. 

END 
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1. Input deck for matrices ^ A^* | an d 

The input of matrices and Jb^*J, a total of 900 elements, 

is accomplished by reading them from a deck of ISO cards, vith 5 
elements on each card, arranged in Format 5F16.7* See Appendix A for 
their proper arrangement. 




Total no. of cards = 180 



2. Identification card 

In columns 1-4, ( NJ 0B3 ) , Punch the integer number of ships to be 

calculated in this particular run, with the 
units digit of the integer in column 4. 
Format l4. 

In columns 5 - 80, (QQ) Punch any identification information desired. 

Example: "SCHUYLER OTIS BLAND TEST RUN." 

Format 19A4. 

Total no. of cards = 1 



3 . Input for Ships 

a. If input Mode 'I is being used for a ship, punch the following 3. cards 
for each ship: 

Card 1 



In columns 1-4, .(INOPT), Punch a positive integer. Example: "+5". 

Format 14. 



Card 2 



In columns 1 - 10, (XLBP), 



Punch the length between perpendiculars in feet. 
Example: "IiOO.O" Format F10.2. 



In columns 11 - 20, (BM), Punch the maximum beam at the load water line in 

feet. Example: "66.0". Format F10.2. 

In columns 21 - 30, (H), Punch the draft in feet. Example: "27-0". 

Format- F10.2. 



In columns 31 - 40, (CB), 



Punch the block coefficient. Example: "0.651". 
Format F10.2. 



In columns 4l - 50, (Vl-EIT), Punch the lowest speed in knots for which you 

desire information on 3HP and 3HP. Example : 
"8.557"- Format F10.2 
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In columns 51 - 60, 



(VKAX), Punch the highest speed in knots for which 
you desire information on EHP and SEP. 

Example : "19.25". Format F10.2. 

In columns 6l - 70, (VIKC), Punch the increment of speed at which you wish 

information on EHP and SEP. Example: "1.07" 
Format F10.2. 



Card 3 

In columns 1 - 10, (DCF), Punch the frictional ship correlation factor 

you desire to use. Example: ".0004". 

Format F10.2. 



In columns 11 - 20, (XLCB), Punch the position of LCB from ^ in feet, 

with a + sign if forward and a - sign if 
aft of ^ . Example: 4.0". Format F10.2. 



In columns 21 - 30, (WSURF), Punch the wetted surface area in square feet. 

Example: "39994 .0". Format F10.2. If wetted 

surface is not known, punch a "0.0" in this 
space and the program will calculate it. 



In columns 3-- ~ 40, (PF°) 



2 , 4 

Punch the density of vater in lb. sec /ft. 
at the temperature desired. Example: "1.9905"* 
Format F10.2. 



In columns 4l - 50, (GNU), Punch the kinematic, viscosity in ft.^/sec. 

Example: "I.2817E-05", for I.2S17 x 10"5. 
Format E10.2. 



Mote : If information on EHP and SEP is desired for only one speed for the 

ship, e.g. service speed, on card 2 punch this speed in the spaces for 
VMM and VMA.X, and punch "0.0" in the space for VINC. 

Total no. of cards = NJOBS x 3- 



b. If input Mode II is being used for a ship, follow the instructions 



for Mode I above with the following exceptions: 



Card 1 

In columns 1-4, (IMOPT), Punch a negative integer. Example: "-5". 

Format 14. 



Card 2 



In columns 11 - 20, (BE), 



Punch the beam to draft ratio. Example: 
"2.5". Format F10.2 



In columns 21 - 30, (DEL), 



Punch the displacement in tons. Example: 
"10000.0". Format F10.2. 
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* 



I 










In columns 31 - 40, (CP), 



Punch the prismatic coefficient. Example 
".612". Format F10.2. 



Card 3 

Identical to Mode I. 

Mode I and Mode II input ships can he mixed at will as long as the 
cards are punched as above. 

Graphs for making corrections to the computer output for variation 
in LCB position are contained in Appendix B. 

4. Loading of cards 

The cards are loaded in the following order : 

1. fein EHP-SH? Program Deck 

2. Subroutine CEF 

3 . Input deck for matrices [A ^ and[B n -j 
k. The identification card 

5. Input cards for ships. 
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